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SECTION  1 
INTRODUCTION 


A  single  area  fire  started  by  a  nuclear  weapon  burst  is  likely  to 
have  a  well-defined  boundary  and  little,  if  any,  increase  in  the  fire 
area  In  the  few  recorded  area  fires,  spread  did  not  occur  [Bond, 

19*16;  Irving,  1963;  and  Miller,  1968];  the  high-velocity  fire  winds 
generated  by  the  strong  buoyancy  of  the  area  fires  [Small  and  Larson, 
1986]  limited  the  outward  spread.  Unless  the  ambient  wind  velocity  is 
greater  than  the  induced  fire-wind  velocity,  the  area  fire  system  is 
likely  to  be  stable  and  not  grow.  There  can,  of  course,  be  exceptions 
due  to  topography,  uneven  fuel  loadings,  or  local  instabilities  at  the 
fire  boundary  (fingering)  that  could  lead  to  changes  in  the 
perimeter.  In  general,  stable  fronts  are  expected  for  single  area 
fires . 

The  fire-wind  velocity  is  greatest  at  the  fire  perimeter  and 
decreases  slowly  with  distance  from  the  fire.  This  is  for  a  mature 
fire— one  that  has  burned  long  enough  to  establish  a  far  field.  (Such 
fires  are  more  likely  to  occur  in  urban  areas  than  in  wildlands  due  to 
the  larger  fuel  concentrations  in  an  urban  area  and  correspondingly 
longer  burning  times.)  Calculations  show  that  velocities  on  the  order 
of  1  m/s  are  possible  at  3  to  *1  fire  radii.  Such  far-field  velocities 
are  large  enough  that  separa*^  area  fires  could  interact  if  spaced 
sufficiently  close  to  each  other. 

In  many  urban  areas,  potential  targets  are  widely  spaoed  (in 
terms  of  weapon  radii)  and  multiple  bursts  are  likely.  Thus,  several 
area  fires  could  develop.  Initially,  each  area  fire  would  burn  inde¬ 
pendently;  however,  as  tho  flow  fields  develop,  interactions  leading 
to  changes  in  the  fire  boundary  are  possible.  The  area  fires  could 
merge;  smaller  or  weaker  fires  migrating  toward  the  larger  or  more 
intense  fires. 

In  this  report,  we  consider  the  interactions  of  multiple  large  area 
fires  separated  by  distances  on  the  order  of  a  fire  diameter,  and 


determine  when  fire  spread  is  probable  for  some  simple  lay  downs.  A 
simple,  analytical  model  of  the  fire  plume  and  its  effects  on  the 
surroundings  is  developed.  Distributions  of  singularities  such  as 
sinks,  sources,  and  vortices  are  then  defined  such  that  the  flow  field 
induced  by  a  single  large  area  fire  [Small,  Remetch,  -•  id  Brode,  1 984] 
is  reproduced.  Although  basically  inviscid,  the  model  effectively 
accounts  for  viscous  effects. 

Next,  two  or  more  such  distributions  are  placed  in  close 
proximity  and  the  interactive  influence  on  the  .'ire  winds  and  the 
spreading  smoke  plume  flow  fields  are  calcula'  ed.  This  simplified 
model  is  used  to  approximate  the  interaction  effects  without  recourse 
to  hydrocode  solution  of  the  three-dimensional  unsteady  Navier-Stokes 
equations.  We  specifically  study  the  ways  in  which  nearby  fires  and 
their  associated  flow  fields  can  promote  the  spreading  of  an  initially 
circular,  uniformly  burning  area.  The  present  results  reproduce,  at 
least  qualitatively,  the  direction  of  changes  in  individual  fire 
boundaries  due  to  neighboring  fires. 
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SECTION  2 

ANALYSIS  OF  AN  ISOLATED  AREA  FIRE 


Large  area  fires,  defined  such  that  the  radius  of  the  burning 
area  is  on  the  order  of  the  atmospheric  thickness  [0(10^  m)],  produce 
a  typical  flow  field  consisting  of  a  ground  level  inflow  of  air,  a 
flame  zone  of  radius  R  and  height  Hf,  a  buoyant  plume  of  combustion 
gases  and  smoke  moving  upwards  (with  possible  crosswind  effects),  air 
entrained  from  around  the  plume,  and  a  spreading  "fountain”  at  the 
plume  equilibrium  height  H  at  which  the  plume  gases  are  in 
equilibrium.  That  pattern  is  obtained  by  relatively  high  resolution 
calculations  of  the  time-dependent,  viscous,  buoyant  flow  field 
[Small,  Remetch,  and  Brode,  1984].  Similar  hydrocode  solutions  have 
shown  relatively  good  agreement  with  observations  and  data  from  scale 
experiments  [Small,  1985]  and  the  Flambeau  tests  [Nielsen,  1970]. 

Ine  flow  field  of  an  isolated  fire  outside  of  the  burning  zone 
(r  >  R)  is  modeled  by  a  vertical  distribution  of  horizontal  planar 
sources  and  sinks  placed  at  the  fire  centerline,  with  sinks  at  the 
lower  parts  to  generate  inflow  and  sources  at  the  top  (see  Fig.  1). 
That  approach  is  consistent  with  the  observation  of  Smith,  Morton,  and 
Leslie  [1975]  and  Cox  and  Chitty  [1980]  of  a  "sink-like"  inflow  at 
ground  level,  for  large  fires.  Swirl  could  be  easily  included  by 
adding  a  vortex  with  vertical  axis.  That  would  not  affect  the  sink 
model  at  ground  level— the  flow  remains  essentially  radial  even  if 
strong  swirl  exists  [Emmons,  1 9651 -  Moreover,  we  have  recently  shown 
that  a  large  plume  is  not  likely  to  swirl  [Weihs  and  Small,  1985]. 

The  singularity  strength  is  modeled  by  a  combination  of  "basic" 
sinks  that  describe  the  fire  requirements  and  "entrainment"  sinks  that 
model  the  added  entrained  (plume)  flow. 

The  strength  of  the  basic  sinks  is  obtained  from  the  energy  per 
unit  time  Qp  produced  by  a  fire  of  radius  R.  This  is 


a.  Axonometric  description  of  fire  and  plume. 
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Typical  streamline 


b.  Vertical  section  through  fire  center  and  description  of 
coordinates,  parameters,  and  model  singularities. 

Note:  Burning  zone  height  is  greatiy  exaggerated  to  show  detail. 


Figure  1.  Schematic  representation  of  large  area  fire. 
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where  J  is  the  heat  addition  rate  per  unit  volume.  If  C  is  the 
average  caloric  value  of  the  combustibles  per  unit  mass,  the  mass  flux 
of  air  required  to  sustain  the  fire  is 


A  =  S 


(2) 


where  S  is  the  air/fuel  ratio. 

Computations  [Small,  Remetch,  and  Brode,  1984]  show  that  the 
inflow  velocity  in  the  area  surrounding  the  fire  (the  outer  flow)  is 
roughly  described  by  an  r"1  dependence,  with  no  noticeable  changes 
over  the  first  100  m  of  altitude.  Thus,  we  assume  the  basic  sinks  to 
be  a  linear  vertical  distribution  (of  an  as  yet  unknown  height)  of 
planar  sinks.  That  implies  each  horizontal  section  is  independent  in 
terms  of  the  horizontal  component  of  flow,  i.e.,  a  source  or  sink  at  a 
given  height  only  drives  the  flow  at  that  height  [Lamb,  1945].  The 
height  h  is  obtained  from  the  computed  ground  velocities  ug  by 
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where  r  is  the  distance  (r  £  R)  at  which  ug  is  measured.  The  basic 
singularity  is  of  strength  B  -  A/h  (in  square  meters  per  second) 
distributed  over  a  vertical  distance  h. 

The  induced  singularities  are  obtained  by  calculating  the  inflow 
or  outflow  at  greater  heights.  To  estimate  the  entrained  flow  we 
first  need  to  model  the  fire  plume.  The  large  fires  described  here 
are  unusual  in  the  sense  that  the  ratio  of  plume  height  to  fire 
diameter  is  0(1)  or  less.  Equilibrium  of  the  plume  is  controlled  by 
crosswinds  or  gradients  in  the  atmosphere,  rather  than  by  entrainment 
of  ambient  air.  Accordingly,  the  point  source  model  [Turner,  1973]  of 


SECTION  3 

INTERACTIONS  BETWEEN  AREA  FIRES 


To  examine  the  possible  spreading  of  area  fires,  we  consider  the 
effect  of  a  neighboring  area  fire  on  the  flow  field  at  the  perimeter. 
Large  area  fires  such  as  those  at  Hamburg,  Dresden,  or  Hiroshima  have 
exhibited  sharply  defined  bounds  of  the  burnt  area,  indicating  that 
strong  inflow  fire  winds  limit  outward  spread.  This  implies  that  for 
some  minimum  radial  local  inflow  velocity,  the  fire  cannot  spread. 

That  minimum  velocity  V0  is  a  function  of  the  energy  density  of  the 
fire  J,  and  possibly  other  parameters  such  as  local  curvature  of  the 
flame  front  1/R,  continuity  of  fuel,  etc.  The  minimum  velocity  V0  is 
a  necessary,  but  not  sufficient  limit  (in  mathematical  terms),  i.e., 
if  0  <  | u j  <  |  V0 1  the  fire  can,  but  does  not  necessarily,  spread. 

For  positive  (u  >  0)  velocities  at  the  perimeter,  we  assume  the  fire 
will  spread. 

As  a  result  of  the  above  arguments,  it  seems  that  the  dominant 
factors  regulating  fire  spread  appear  at  ground  level,  within  a  layer 
of  thickness  similar  to  that  of  the  fire  region.  As  shown  previously, 
the  flow  in  the  bottom  layer  is  essentially  that  of  a  two-dimensional 
sink.  That  will  now  be  applied  to  the  interaction  problem.  Looking 
first  at  fires  that  are  established  simultaneously,  we  consider  two, 
nonequal  size  fires  of  sink  strengths  B]  and  B2,  respectively,  where 
B1  *  B2. 

For  simplicity,  we  define  the  length  unit  as  the  larger  fire 
radius  (i.e.,  R-]  =  1).  Then  R2  =  b  S  1 ;  the  equality  sign  stands  for 
equal-sized  circular  fires.  The  velocities  induced  by  each  of  the 
fires  independently  follow  the  form 


il 

2irr 
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(ri  *  Ri) 


(5) 


where  r^  represents  the  radial  distance  from  the  respective  fire 
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center.  Thus,  the  speed  at  a  point  P  on  the  perimeter  of  the  smaller 
fire  on  the  line  connecting  the  centers  (see  Fig.  3)  is 


p  2ir(D  -  b)  2trb  ’  v  ; 

where  D  is  the  distance  between  fire  centers,  normalized  by  R.  When 
Up  £  0  the  smaller  fire  tends  to  spread  toward  the  larger  fire.  That 
process  is  self-amplifying.  The  outward  velocity  of  the  moving  area 
fire  is  larger  because  it  is  closer  to  the  "dominant"  fire.  For  the 
limiting  case  (up  =  0), 


1  D  -  b  D 

—  =1  1  33  — 

B2  b  b 


D  »  b  1  +  , 

m  B2 


where 


B~  *  15  Dm  *  1  * 


When  the  fires  are  initially  of  equal  strength  (B-|  =  B2)  Dm  =>  2,  the 
fires  merge  only  if  initially  touching.  For  all  other  cases,  there  is 
a  finite  distance  over  which  the  fires  tend  to  combine.  The  result 
depends  on  the  fire-induced  winds  only.  Ambient  winds  may,  of  course, 
affect  the  possibility  of  outward  spread. 

A  rough  estimate  Can  be  determined  of  the  distance  Dm,  as  a 
function  of  the  ratio  of  fire  sizes,  we  assume  both  fires  to  burn 
similar  terrain  or  urban  area,  so  that  the  energy  addition  J  and  flame 
height  Hf  are  equal  for  both.  Then,  the  sink  strength  B  from  Eqs.  (1) 
and  (3)  is  proportional  to  the  radius  squared,  i.e., 
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(10) 


or  from  Eq.  (9) , 


b  +b 


(11) 


The  limit  b— »0  implies  that  a  point  fire  would  spread  an  infinite 
distance  since  the  large  fire  serves  as  an  external  wind  driving  the 
small  fire.  For  the  quasi-steady  model  to  be  relevant,  the  speed 
induced  at  point  P  (Fig.  3)  should  at  least  be  larger  than  the  local 
turbulent  fluctuations  (ambient  wind  or  fire  generated).  Such 
velocity  fluctuations  are  a  higher  fraction  of  the  inflow  speed,  the 
smaller  the  fire.  Thus,  we  assume  that  b  should  be  greater  than  0.1. 

The  changes  in  area  burned  due  to  multiple  fire  interactions  can 
be  obtained  by  comparing  areas  for  the  two  fires  separately  (tt  +  irb2) 
and  the  area  included  between  them  (see  Fig.  3a).  The  ratio  is 


F  .  . .  a  -  +  -  - 

1  .  .2,  2  ir  4  .  u2  ’ 


(12) 


ir(  1  +  b  ) 


1  +  b 


and,  using  Eq.  (11) 


1  +  b  +  1 

2  irb 


(13) 


That  can  be  seen  as  an  upper  bound  for  the  interaction,  accurate  when 
b  — >1  • 

When  the  ratio  of  initial  fire  radii  is  small,  a  different  flow 
pattern  can  develop.  The  Inflow  toward  the  larger  fire  is  radial  with 
respect  to  its  center  and  when  the  secondary  fire  is  much  smaller,  a 
narrowing  spread  area  develops  (Fig.  3b).  That  area  is  approximately 
bD  -  ir  (tan-1  b/D)  (when  b  <’  0.5)  so  that  the  ratio  of  the  merging 
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system  to  the  initial  area  fires  is 


bD  +  ir(l  -  tan 

F  *  — - L__ - U  .  (14) 

ir(1  +  b  ) 

The  approximatibn  results  from  neglecting  parts  of  the  smaller  fire 
circle.  The  approximation  is  best  for  small  b,  and  in  Fig.  4,  F2 
describes  the  limit  for  b  «  1.  Substituting  Eq.  (11)  into  Eq.  (14) 
we  obtain 


(15) 


where  F2  tends  to  a  finite  limit  (»  1.32)  for  b  =  0  (see  Fig.  4).  The 
actual  area  ratio  will  probably  fall  between  the  curves  F \  (for  b  «  1) 
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Note:  Curve  Fi  describes  pattern  of  Fig.  3a,  and  curve  F2  shows  pattern  of  Fig.  3b. 
The  solid  line  interpolates  between  limiting  cases. 

Figure  4.  Ratio  of  total  system  fire  area  to  sum  of  initial  area 
fires  F,  as  function  of  ratio  of  individual  fire  radii. 


13 


i 


and  F2  for  (b — -0).  Thus,  an  increase  of  up  to  50  percent  in  area 
burned  due  to  the  interaction  is  possible  for  fires  with  a  size  ratio 
of  approximately  3,  separated  by  =  3.5  large  fire  radii  (shown 
schematically  in  Fig.  3a-b) . 

The  results  shown  in  Fig.  4  consider  the  merging  of  area  fires 
when  a  net  nonzero  velocity  is  established  between  than.  This  is 
spread  in  the  wind  direction  by  the  weaker  fire;  one  area  fire  is 
stationary.  Fires  can  also  spread  against  the  wind  when  the  fire- 
induced  speed  is  less  than  some  minimum  V0.  That  velocity  is  a  func¬ 
tion  of  the  fire  shape  and  strength.  The  strength  influences  the 
updraft  and  entrainment  at  the  perimeter,  which  in  turn  can  sweep  up 
fire  brands  and  drop  them  beyond  the  fire  area.  The  mechanisms  deter¬ 
mining  a  minimum  value  for  V0  are  not  well  established.  We  thus  vary 
V0  parametrically  and  examine  the  area  fire  interactions.  The  inflow 
velocity  U2  at  the  perimeter  of  an  isolated  fire  of  radius  R2  a  bRi  is 
assumed  related  to  V0  by  Ku2-  Therefore,  from  Eq.  (5) 


Fran  Eq.  (6),  with  V0  =»  up,  we  obtain,  after  some  manipulation, 


1  +  (K  +  1 )b 

Dm  ■  '  °7' 

which  reduces  to  Eq.  (11)  when  K  =  0,  as  expected.  Essentially,  a 
finite  value  of  V0  does  not  result  in  any  qualitative  changes  The 
range  over  which  interactions  are  possible  is  increased  and  similarly, 
the  spread  area  is  increased  (see  Fig.  5). 

Since  area  fires  have  a  life  span  limited  by  the  amount  of  com¬ 
bustibles  available,  there  is  a  different  limitation  on  the  increase 
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Note-.  The  parameter  K  is  ratio  of  incoming  velocity  a  fire  can  overcome  over  inflow 
velocity  at  fire  perimeter.  Straight  tine  indicates  distance  Dm  for  initially 
touching  fires. 


Figure  5.  Changes  In  maximum  nondimensional  distance  between  fire 
centers  as  function  of  fire  radii  ratio. 
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of  burn  area.  The  larger  fire,  assumed  fixed  in  size*,  will  have  an 
approximate  lifetime  of 


T 


1 


(18) 


which  is  the  time  required  to  turn  all  the  combustible  material  in  the 
fire  area. 

The  spread  rate  of  the  fire  at  any  location  is  at  most  ap¬ 
proximately  equal  to  the  airspeed  at  that  location.  On  the  line 
connecting  the  fire  centers,  which  is  the  line  of  fastest  spread,  the 
speed  is,  from  Eq.  (5) 
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2irr 


B 


2 


2ir(D  -  r)  ’ 


(19) 


For  a  rough  estimate  of  the  time  required  for  the  fires  to  merge,  one 
integrates  Eq.  (19)  as 


(20) 


where  r  -  r(t)  such  that  r(0)  -  D  -  b.  The  result  is  only  approximate 
sinoe  the  shape  of  the  merging  fires — especially  that  of  the  smaller 
fire--is  no  longer  circular.  The  induced  velocities  can  no  longer  be 
modeled  by  simple  sinks.  A  numerical  marching  procedure  that  adjusts 
the  sink  distribution  to  acoount  for  ohanging  (e.g.,  elliptical)  fire 
areas  can  be  used.  When  i  i  D  -  b  -  1 ,  where  the  large  fire  initial 

^Actually,  once  the  initially  smaller  fire  has  grown,  ar.d  if 
the  initial  ratio  b  was  close  to  unity,  it  might  become  the  dominant 
fire  and  attract  the  other.  This  does  not  change  the  results  shown 
in  Figs,  t  and  5. 
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radius  is  unity,  the  results  illustrated  in  Figs.  H  and  5  are 
obtained.  Otherwise,  only  partial  covering  of  the  areas  shown  in 
Fig.  3  is  achieved.  As  the  initial  areas  burn  out,  the  intermediate 
area  continues  to  spread. 

In  the  case  of  more  than  two  fires  burning  simultaneously,  the 
point  or  points  at  which  the  velocity  on  the  perimeter  changes  direc¬ 
tion  is  not  known  a  priori,  as  it  depends  on  the  two-dimensional 
positions  of  the  flames,  as  well  as  their  strength.  The  velocity  at 
any  point  can  be  estimated  by  considering  the  following  superposition 
of  singularities 


B2  B3 

- - -  f  - 


2ur1  2ir(r  -  S2)  2ir(r  -  S  ) 


+  •  .  *  + 


2ir(r  -  S  ) 
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i=i 


V  ’ 


(21) 


where  r  is  the  distance  from  the  point  in  question  to  the  center  of 
the  fire  designated  as  the  "first"  fire.  The  choice  is  arbitrary  and 
can  be  made  by  size  or  geographical  location. 
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SECTION  4 
RESULTS 


In  this  section,  we  apply  the  singularity  theory  to  several  com¬ 
binations  of  area  fires  and  calculate  the  interactions.  A  10-km 
radius  fire  serves  as  a  basis  for  obtaining  parameters  for  a  single 
area  fire  first,  and  then  for  determining  multiple  area  fire  interac¬ 
tions  (see  Sec.  2) . 

We  first  consider  whether  the  flow  field  defined  by  the  calcula¬ 
tions  can  be  assumed  quasi-steady.  Computed  ground  level  flows  are 
10  m/s  at  the  fire  perimeter  (Fig.  2)  and  decay  roughly  as  r-1  with 
distance  from  the  fire  center.  Thus,  far-field  inflow  speeds 
of  0(1  m/s)  are  expected,  and  for  distances  of  0(10^  m),  time  scales 
of  10^  s  are  obtained. 

The  fire  [calculated  by  Small,  Remetch,  and  Brode  (1984)]  i3 
assumed  to  burn  in  an  area  with  fuel  loading  between  that  of  the 
average  urban- suburban  combustible  load  of  40  kg/m2  [National  Academy 
of  Sciences,  1985]  and  a  city  such  as  Hamburg  for  which  the  load  was 
estimated  at  157  kg/m2  [Carrier,  Fendell,  and  Feldman,  1985].  The 
caloric  equivalent  of  those  combustibles  is  1.86  •  10?  J/kg  so  that 
the  total  heat  released  from  the  fire  per  square  meter  is  between 
7,44  •  10®  and  29.2  *  10®  J/m2.  The  averaged  heat  addition  rate  was 
1  kW/m3  over  a  height  of  100  ra  so  that  each  square  meter  of  surface 
area  accounts  for  100  kW/m2,  Thus,  the  fire  duration  is  at  most  from 
0.744  to  2.92  •  10^  s,  i.e.,  between  2.07  and  8.1  h,  depending  on  the 
type  of  area  burned.  The  fire  duration  is  of  the  same  order  of  mag¬ 
nitude  as  the  characteristic  times  so  that  a  steady  flow  is  not  fully 
established.  Nevertheless,  "persistent”  (i.e.,  qua3i-steady ,  slowly 
changing)  patterns  appear  after  -  3000  s.  Also,  the  large  size  of  the 
aggregate  fire,  which  is  actually  composed  of  a  multitude  of  smaller 
fires  is  assumed  to  cause  averaging  of  locally  periodic  buoyant  ef¬ 
fects  [Turner,  1983].  As  a  result,  the  velocity  field  exterior  to  the 
fire  behaves  approximately  like  a  steady  flow. 


Figure  2  shows  the  radial  inflow  calculated  at  60  and  90  min 
after  fire  initiation.  To  determine  the  basic  sink  distribution 
strength,  we  use  Eqs.  (1)  through  (3),  to  calculate  the  mass  and  thus 
the  sink  strength.  For  R  =  10^  m,  Hf  =  102  m,  J  =  1  kW/m3,  and 
C  -  1.86  •  10?  J/kg,  Q  =  3.1  ^  •  10^  j/s  and  A  =  2.53  •  10?  kg/s. 
Taking  the  values  of  the  radial  inflow  at  3600  and  4800  3  after  the 
fire  start,  and  an  incoming  air  density  of  p  =  1 .2  kg/m3,  the  equiv¬ 
alent  sink  height  is  h  *  30  m.  That  value  describes  the  region  where 
u  is  proportional  to  r-1 ,  and  separates  the  effects  of  the  basic  sink 
distribution  from  those  of  the  secondary  vortices  (which  appear  as  a 
"hump"  in  the  velocity  distribution  shown  in  Fig.  2).  This  hump  is 
the  footprint  of  a  concentrated  vortex  ring  that  is  moving  radially 
outward  as  time  progresses.  The  biggest  increment  is  observed 
directly  below  the  vortex  ring  "core,"  falling  off  as  r-1  cos  0,  where 
0  is  the  angle  between  the  vertical  and  the  line  connecting  each  point 
and  the  vortex  core.  In  keeping  with  our  model  of  superposed  sin¬ 
gularities,  the  vortex  strength  is  calculated  separately  below. 

The  sink  height  (h  =>  30  m)  necessary  to  provide  the  air  required 
for  the  fire  is  three  times  less  than  the  100  m  burning  zone  height, 
and  within  the  zone  over  which  the  inflow  velocity  does  not  change 
with  height.  The  result  is  consistent  with  the  assumption  of  a 
uniform  basic  sink  strength.  The  sink  strength  B  per  unit  height,  is 
B  =  8.43  •  i05  m2/s.  Thus,  the  fire  has  more  than  300  percent  excess 
air,  so  that  burning  can  be  unhampered  and  complete  at  a  height  of 
100  m. 

An  alternate  estimate  of  the  sink  strength  can  be  obtained  by 


kinematic  means  [Eq.  (3)],  directly  from  the  velocity  curves 
(outside  the  fire)  of  Fig.  2.  For  an  elapsed  time  of  36OO  s, 

B  -  6.94  •  10^  m2/s  (i.e.,  within  18  percent),  and  at  t  -  4500  s, 

B  =  5.26  •  105  m2/s.  Some  of  the  discrepancy  might  result  from  our 
neglect  of  unsteady  effects  such  as  growth  and  movements  of  annular 
vortices.  For  times  shorter  than  1  h,  however,  the  starting  transient 
dominates.  External  speeds  of  <  1  5  m/s  were  calculated  so  that  at  1  h 
the  effect  of  the  fire  has  only  reached  <5  •  10**  m.  This  can  be 
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clearly  seen  in  the  plot  of  versus  r  for  t  =  2700  where  the  effects 
of  the  fire  have  only  reached  =  4  •  101*  m  (see  Fig.  6).  In  any  case, 
agreement  between  the  two  independent  estimates  is  encouraging. 

The  flow  around  the  basic  sink  if)  now  separated  into  induced 
sinks  and  a  vortex.  The  sink  distribution  above  the  fire  is  again 
obtained  by  plotting  the  horizontal  inflow,  this  time  at  various 
heights.  For  the  range  0  to  3000  m,  the  horizontal  velocity  exhibits 
an  r-1  dependence  on  r  >  10^  m.  The  results  are  given  in  Fig.  7. 

The  vortex  models  the  turning  of  lower-level  inflow  into  upper- 
layer  spreading.  It  is  essentially  a  vortex  ring  of  radius  =■  18  km  at 
60  min  and  =  36  km  at  90  min,  at  a  height  of  =■  3  km.  The  vortex 
height  is  approximately  halfway  between  the  ground  and  the  major 
outflow  layer,  with  additional  less  well-defined  vortex-like  struc¬ 
tures  above  it.  Thus,  the  vortex  ring  can  be  modeled  as  being  con¬ 
strained  between  two  horizontal  planes.  This,  in  turn,  is  mathemati¬ 
cally  equivalent  to  having  an  infinite  row  of  image  vortex  rings  of 
alternating  direction  of  rotation  with  a  common  axis,  equally  spaced 
vertically  at  about  6  km.  The  vortices  included  in  the  infinite  row 
are  essentially  reflections  of  the  original  ring  in  the  two  constrain¬ 
ing  plumes.  The  velocity  induced  by  the  vortex  on  the  ground  is  used 
to  establish  the  vortex  strength  (see  Appendix  A). 

The  vortex  strength  at  3600  s  is  r  =  1600  m2/s  and  its  contribu¬ 
tion  to  the  velocity  field  is,  at  most  about  12  percent  of  the  ground 
velocity.  Figure  8  illustrates  the  ground-level  inflow  velocity 
variations  from  the  sink-flow  pattern,  clearly  pinpointing  both  the 
vortex  (the  peak)  and  the  fire  edge  effect.  The  outer  flow  is  incom¬ 
pressible,  so  that  the  increment  in  ground-level  velocity  results  from 
mass  coming  from  above  while  moving  in  and  moving  back  upwards  at 
radial  distances  less  than  the  instantaneous  vortex  position.  The 
difference  between  the  calculated  ground  flow  and  sink  flow  at  4800  s 
is  that  the  vortex  has  moved  out  to  -  36  km  and  another  vortex  is 
forming  at  15  to  16  km  from  the  fire  center  (still  relatively  weak). 
The  strength  of  the  original  vortex  is  unchanged  (its  height  is  con¬ 
stant  and  the  maximal  induced  velocity  has  changed  by  less  than 
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5  percent  to  slightly  above  0.8  m/s ) .  The  center  of  vorticity  is  mo  zing 
radially  outward  at  an  average  speed  of  10  m/s  during  the  period  3600 
to  4800  s  after  inception  (an  indication  of  the  growth  of  the  cloud  of 
combustion  products  near  the  tropopause). 

Figure  9a  illustrates  the  ground-level  flow  field  generated  by  a 
configuration  of  three  fires  all  burning  an  area  of  uniform  combus¬ 
tible  density.  The  strength  (mass  burned  per  unit  time)  of  the  two 
smaller  fires  (Nos.  2  and  3)  is  36  percent  of  the  larger  fire  (propor¬ 
tional  to  the  radii  ratio  squared).  The  streamlines  between  fires  1 
and  2,  and  1  and  3,  respectively,  are  similar,  and  suggest  a  burn 
pattern  such  as  sketched  in  Fig.  3.  In  contrast,  Fig.  9b  illustrates 
a  case  in  which  fire  2  is  attracted  to  fire  1 ,  but  fire  3  is  not 
attracted  to  fire  1.  Thus  a  burning  corridor  is  established  only 
between  fires  1  and  2.  That  illustrates  the  effect  of  intensity  on 
merging  of  area  fires. 

Figures  9c-d  show  a  different  placement  of  the  three  bursts  (and 
fires)  described  in  Figs.  9a-b.  In  Fig.  9d,  fire  3  burns  a  similar 
size  area  but  with  twice  the  combustible  loading.  Fire  2  is  attracted 
to  fire  3  and  not  to  the  larger  fire  1.  Those  results  illustrate  the 
influence  of  parameters  such  as  fire  areas,  relative  positions,  and 
fire  strengths  on  the  stability  of  an  area  fire  or  system  of  area 
fires.  The  configuration  shown  in  Fig.  9a  leads  to  a  superfire  encom¬ 
passing  all  the  initial  fires.  Alternatively,  three  fires  (Fig.  9c) 
similar  in  size  to  those  in  Fig.  9a,  remain  separate.  In  those  cases, 
the  relative  positions  determine  whether  the  fires  merge.  The  com¬ 
parisons  of  both  Figs.  9a-b  and  Figs.  9c-d  illustrate  the  influence  of 
area  fire  intensity  or  strength.  A  heavier  fuel  loading  in  one  dis¬ 
trict  may  cause  merging.  Similarly,  different  burst  times  may  lead  to 
merging  due  to  unequal  (at  least  initially)  burn  rates.  Those  results 
suggest  that  multiburst  optimization  procedures  can  be  applied  to 
maximize  or  minimize  the  burned  area. 

It  should  be  mentioned  that  the  calculations  shown  in  Figs.  9a-d 
are  instantaneous  pictures  of  each  fire  flow  field.  After  two  or 
more  fires  have  produced  burning  corridors  or  merged,  the  total  sink 


Burst 

X-center 

Y^center 

Radius 

Strength 

1 

0.00 

0.00 

1.00 

1.00 

2 

2.00 

0.00 

0.58 

0.33 

3 

•1.50 

•1.50 

0.58 

0.67 

b.  Fires  t  and  2  merge;  fire  3  remains  separate. 
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Figure  9.  Ground  level  flow-field  streamlines  for  three  fires  in 
proximity  (Continued). 
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Note:  Fires  2  and  3  have  radii  60  percent  that  of  fire  1,  and  all  burn  in  a 
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Figure  9.  Ground  level  flow-field  streamlines  for  three  fires  in 
proximity  (Continued). 
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d.  Fires  2  and  3  merge;  fire  1  remains  separate. 

Note:  Fire  3  now  burns  twice  the  combustible  per  unit  area  and  time 
as  the  other  fires. 

'Figure  9.  Ground  level  flow-field  streamlines  for  three  fires  in 
proximity  (Concluded). 


strength  changes,  which  can  result  in  a  change  in  the  flow  (stream¬ 
line)  pattern  and  thus  the  burning  pattern.  An  integration  procedure 
can  be  devised  to  account  for  the  quasi-steady  changes  in  the  fire 
sizes  and  shapes. 

All  the  above  calculations  assume  V0  =  0.  A  rough  estimate  of 
the  maximum  inflow  velocity  V0  which  a  fire  can  overcome  and  spread 
upwind,  can  be  drawn  by  using  the  criteria  suggested  in  Glasstone  and 
Dolan  [1977].  The  best  current  estimate  requires  a  wind  of  less  than 
8  mph.  That  translates  to  less  than  3.6  m/s  or  K  <  0.28,  i.e.,  within 
the  parameter  range  given  in  Fig.  5. 
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SECTION  5 
CONCLUSIONS 


Hydrodynamic  considerations  show  that  large  area  fires  can  inter¬ 
act  when  in  proximity.  Single  large  area  fires  are  usually  charac¬ 
terized  by  a  fixed  position  and  a  quasi-steady  radial  inflow  of 
ground-level  air.  The  flow  field  inhibits  spreading  in  a  solitary 
fire,  but  when  two  or  more  unequal  fires  are  within  a  few  radii  of 

u.ch  other,  the  larger  or  more  intense  fire  can  "attract"  smaller  or 
weaker  fires,  causing  (for  two  fires)  an  increase  of  up  to  50  percent 
in  total  burned  area.  That  effect  is  greatest  when  the  ratio  of  fire 
sizes  is  about  3,  or  when  they  are  -  3  large  fire  radii  apart.  That 
interaction  can  be  multiplied  if  more  than  one  small  fire  is  near  a 
large  fire.  Similar  interactions  occur  for  fires  of  equal  size,  but 
with  different  burning  ratios  (combustible  loads). 

Our  conclusions  hold  for  simultaneously  burning  fires.  Different 
merge  criteria  apply  for  fires  starting  to  burn  in  the  vicinity  of  an 
established  large  area  fire.  Preliminary  considerations  suggest  that 
the  configurations  similar  to  those  illustrated  in  Fig.  3b  will  be 
obtained,  with  much  larger  separations  still  leading  to  fire 
coalescence.  Such  time-dependent  effects,  as  well  as  area  fire  con¬ 
figurations  representative  of  lay  downs  in  actual  urban  areas  will  be 
considered  in  a  future  PSR  study. 
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APPENDIX  A 

CALCULATION  OF  VORTEX  STRENGTH 


The  strength  of  the  vortex  used  to  model  the  flow  outside  of  the 
fire  plume  is  obtained  by  an  inverse  technique  using  the  velocity 
induced  at  the  ground.  The  stream  function  describing  flow  outside  of 
an  isolate'’'  vortex  ring  is  [Lamb,  1932] 

*  -  ^  (r!  +  »'2)[F1  (X)  -  E 1  ( X )  ]  ,  (22) 


where  r-|  and  rg  are  the  least  and  greatest  distances  frcm  the  point  in 
question  to  elements  of  the  ring,  X  =  ^2  ~  ri  (r2  +  r^ ) ,  and  F-|  and  E-| 
are  the  complete  elliptic  integrals  of  the  first  and  second  kind, 
respectively.  Equation  (22)  may  now  be  used  as  a  basis  for  generating 
the  velocity  due  to  the  infinite  column  of  image  vortices  by 
summation.  The  radial  velocity  at  the  ground  is  obtained  as  the 
derivative 
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1  <N> 
r  5)z 


(23) 


A  great  simplification  in  the  calculation  is  made  by  assuming 
that  the  flow  field  directly  beneath  the  ring  vortex  is  a  rectilinear 
vortex.  That  assumption  introduces  only  small  inaccuracies,  as  the 
radius  of  the  ring  vortex  (-  15  to  30  km)  is  much  greater  than  the 
vertical  distance  of  the  point  in  question  -  3  km  (directly  beneath 
one  element)  and  the  error  is  proportional  to  the  inverse  of  the 
r.itio.  The  simplified  calculation  is  presented  here  to  describe  the 
process  used  for  the  ring  vortices. 
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The  complex  potential  of  a  single  line  vortex  at  point  (0,  ia)  is 
-irHn(z  -  ia)  where  z  is  the  complex  variable  z  =  x  +  iy.  The  images 
appear  as  vortices  of  opposite  sign  at  (0,  -  ia  +  4  na)  and  of  the 
same  sign  at  (0,  ia  +  4  na)  where  n  goes  from  -®  to  +«  so  that  the 
complex  potential  is 

oo 

W  =  -  [i r  in(z  -  ia  -  4na)  -  ir  £n(z  +  ia  -  4na)]  , 
n=-® 

.  „  „  ,  ,  irz  „  ,  ,  ir(z  -  2ia) 

=■  -ir  in  sin  h  r-  -  in  sin  h  — — r - -  , 

4a  4a 


-if  in 


sin  h  r- 
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sin  h 
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(24) 


The  velocity  induced  directly  beneath  the  vortex  center  is 
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(25) 


The  induced  velocity  is  57  percent  higher  than  for  a  single  vortex. 
That  factor  can  bo  alternatively  understood  as  resulting  from  "squeez¬ 
ing"  the  flow  by  constraining  it  to  move  parallel  to  the  ground.  The 
greatest  velocity  increment  is  about  0.84  m/s  (Fig.  8)  so  that  for  the 
height  of  a  3000  m,  we  obtain  r  ■  1600  m2/s  [Eq.  (25)]. 


34 


APPENDIX  B 
LIST  OF  SYMBOLS 


A  =  air  mass  inflow  rate  [MT-2] 

B  =  sink  strength  per  unit  length  [M/LT] 
b  =  ratio  of  radii  of  small  fire  to  large  fire 
C  =  average  caloric  value  of  combustibles  [L2T~2] 

D  =  maximum  distance  for  definite  spreading  [L] 

Dm  =  maximum  distance  for  definite  spreading  [L] 

F  =  defined  in  Eq.  (13) 

h  =  height  of  basic  sink  distribution,,  simulating  air 
inflow  [L] 

K  =  plume  final  height  [L] 

Hf  =  flame  zone  height  [L] 

J  =  heat  addition  ratio  per  unit  time  and  volume  [ML~^T~3] 

K  =■  defined  by  Eq.  (15) 

m  =  combustible  mass  per  unit  area  burned  [ML"2] 

Qf  =>  total  energy  converted  per  unit  time,  by  fire.  [ML2T~3] 
r  =  radial  coordinate  [L] 

R  =  flame  zone  radius  [L] 

S  =  air/fuel  ratio  of  fire  [0] 
t  =>  time 

u  =  radial  speed 

Ug  =  radial  speed  at  ground  level  [L/T] 

Up  =  velocity  at  point  p  Csee  Eq.  (6)]  [L/T] 

Vf  -  flame  zone  volume  [l3] 
z  =  vertical  coordinate  [L] 

0  »  angular  coordinate 
t  =  fire  lifetime 
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